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It is a well-known fact that, under mild sampling conditions, the restricted Delaunay trian- 
gulation provides good topological approximations of 1- and 2-manifolds. We show that this is 
not the case for higher-dimensional manifolds, even under stronger sampling conditions. Specifi- 
cally, it is not true that, for any compact closed submanifold M of K™, and any sufficiently dense 
uniform sampling L of M, the Delaunay triangulation of L restricted to M is homeomorphic 
' to M, or even homotopy equivalent to it. Besides, it is not true either that, for any sufficiently 

, dense set W of witnesses, the witness complex of L relative to W contains or is contained in 

^ ' the restricted Delaunay triangulation of L. 

Cj ' 

1 Background and definitions 

> , 

\& \ All manifolds considered in this paper are compact closed submanifolds of Euclidean space M n . 

■ The reach of a manifold M, or rch(M) for short, is the minimum distance of a point on M to the 

medial axis of M. All our manifolds have a positive reach. This is equivalent to saying that they 
are ^-continuous, and that their normal vector field satisfies a Lipschitz condition. 
' Given a (finite or infinite) subset L of a manifold M, and a positive parameter e, L is an 

e-sample of M if every point of M is at Euclidean distance at most e to L. In addition, L is 
e-sparse if the pairwise Euclidean distances between the points of L are at least e. Note that an 
e-sparse sample of a compact set is always finite. Parameter e is sometimes made adaptative in 
^ ■ the literature pQ, its value depending on the distance to the medial axis of the manifold. In this 

context, our e-samples are called uniform e-samples. 

For any finite set of points L C ffi n , D(L) denotes the n-dimensional Delaunay triangulation 
of L, and T> M (L) its Delaunay triangulation restricted to a given subset M of M. n . By definition, 
T> M {V) is the nerve of the restriction of the Voronoi diagram of L to X. For any simplex a of T>(L), 
V(cr) stands for the face of the Voronoi diagram of L that is dual to a. The following result comes 
from P|2]: 

Theorem 1.1 If M is a smooth curve in the plane or a smooth surface in 3-space, and if L is a 
finite e-sample of M, with e < 0.1 rch(M), then T> M (L) is homeomorphic to M. 

Let L, W be two subsets of IR n , such that L is finite. Given a point w G W and a simplex 
a = [Po> ' " ' iVk] with vertices in L, w is a witness of o (or simply w witnesses a) if po, ■ • • ,pk 
are among the k + 1 nearest neighbors of w in the Euclidean metric, that is: \/p £ {po, ■ ■ ■ ,Pk}, 
Vg £ L \ {po, ■ ■ ■ ,Pk}, \\w — p\\ < \\w — q\\. The witness complex of L relative to W, or C W (L) for 
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Figure 1: Left: 2-d version of hypersurface M (bold), denned as the boundary of the Minkowski 
sum of hypercube [— -j? (solid) with the ball of radius — centered at the origin (copies of this 
ball are dashed). Hypercube [—A, A] 2 is marked by dotted lines. Right: M and its medial axis. 

short, is the maximum abstract simplicial complex with vertices in L, whose faces are witnessed 
by points of W. From now on, W will be referred to as the set of witnesses, and L as the set of 
landmarks. As pointed out in 0[H], when W samples a manifold M, the witness complex C w (L) 
can be viewed as a discrete version of the restricted Delaunay triangulation T> M (L), and as such it 
should be closely related to it. This is true indeed for curves and surfaces, as stated in the following 
result of 0152]: 

Theorem 1.2 There exists a positive constant c such that, if M is a smooth curve in the plane 
or a smooth surface in 3-space, and if L is an e-sample of M, with e < c rch(M), then C w (L) is 
included in T> M (L) for any set of witnesses W C M, and it coincides with T> M (L) if'W = M. 

2 Negative results 

In this section, we prove that Theorems 11.11 and 11.21 do not hold as is for higher-dimensional 
manifolds. We first show that T> M {L) is not always homeomorphic to M, even though L is an f2(e)- 
sparse 0(e)-sample of M, for arbitrarily small e (Theorem 12. lj) . Our proof builds on an example 
of [6l §11], which deals with hypersurfaces in IR 4 . The intuitive idea is that, when T> M (L) contains 
badly-shaped tetrahedra, called slivers, it is possible to make its normals turn by a large angle 
(say ?) by perturbing the points of L infinitesimally. It follows that the combinatorial structure 
of T> (L) can be modified by small perturbations of M. We then extend our counter-example to 
show that V M (L) may even not be homotopy equivalent to M (Theorem 12. 2p . Finally, we show 
that C W (L) may not be included in V M (L), even for arbitrarily dense sets W C M (Theorem 12.31) . 
The fact that C w (L) may not contain T> M (L) if W C M has already been proved in [9]. 

Theorem 2.1 For any positive constant jj, < | ; there exist a compact closed hypersurface M in M 4 
and an fl(e) -sparse 0(e) -sample L of M, with e = \i rch(M), such that T> M (L) is not homeomorphic 
to M . The constants hidden in the Q, and O notations are absolute and do not depend on \i. 
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Figure 2: Left: tetrahedron [u,v,w,pq] and its dual Voronoi edge. Right: after perturbation of M. 

Proof. Let A = jjj. In R 4 , endowed with an orthonormal frame (x,y,z,t), we construct a 
hypersurface M of reach ^ = -■ Consider the Minkowski sum of hypercube [— -y, — ] 4 with the ball 
of radius ^ centered at the origin. The result is a smoothed-out version of hypercube [—A, A] 4 , as 
illustrated in Figure Q] (left). Let M be its boundary. The reach of M shown in Figure CD 

(right). Let e = [/, rch(Af) = 1, and let d > be an arbitrarily small parameter. Consider points 
u = (1,0,0, A), v = (1,1,0, A), w = (0,1,0, A), and p = (0,0, 5, A). Let c = {\,\, |,A). It is 
easily seen that cq is the circumcenter of [it, v, w, po]. Moreover, all these points belong to M, which 
coincides with hyperplane t = A in their vicinity. Let vq = \\co—u\\ = \\cq— v\\ = \\cq— w\\ = \\cq — po\\. 
We generate an e-sparse 2e-sample Lq of M by an iterative process, starting with Lq = {u, v, w,po}, 
and inserting at each iteration the point of M lying furthest away from the current point set Lq, 
until the farthest point of M is no farther than 2e from Lq. Since M is compact, the process 
terminates, and the outcome is a 2e-sample of M. Moreover, since u,v,w,po lie at least e away 
from one another, and since every point inserted in Lq lies at least 2e away from Lq at the time 
of its insertion, Lq is e-sparse. Finally, no point of ball B(cq,tq) lies farther from {u, v, w,po} 

than 2tq = 2\J | + which is less that 2e since 5 is arbitrarily small. It follows that the interior 

of B(co,ro) contains no point of Lq, which implies that [u,v,w,po\ belongs to V m (Lq), its dual 
Voronoi edge intersecting M at cq. Observe also that, since u,v,w,pq belong to hyperplane t = A, 
the normal of [it, v, w,pq] is aligned with vector (0, 0, 0, 1), as shown in Figure [2] (left). 

We now deform M slightly and create a small bump at cq, such that the top of the bump is 
moved by | into the t-dimension, outward the hypercube. This bump changes the local feature 
size of M. However, since 5 is arbitrarily small, the radius of curvature of the bump can be forced 
to be at least ^, which implies that the reach of M remains equal to — = ^. Furthermore, since 
Co is the center of a Delaunay ball of radius greater than 1, we can assume that 5 is small enough 
for the points of Lq to remain on M. Let c = (^, ^, |, A + |) be at the top of the bump. Since 
the points of Lq are located in hyperplane t = A in the vicinity of [u,v,w,pq], c is equidistant to 
u, v, w,pq, and closer to these points than to any other point of Lq. This implies that the open ball 
B c = B{c, ||c — u||) contains no point of Lq and has u, v, w,pq on its bounding sphere. Hence, B c is 
a Delaunay ball circumscribing [u,v,w,pq], and c belongs to the Voronoi edge dual to [u, v, w,po\. 
Moreover, since u,v,w and (0,0,0, A) are cocircular, dB c passes also through (0,0,0, A). 
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We deform M further by creating another small bump, at point (0, 0, 0, A) this time, so as to 
move this point by 5 into the ^-dimension, outward the hypercube. Let p = (0, 0, 0, A + 6) be the 
top of the bump — see Figure [2] (right). A quick computation shows that ||c — p\\ = ||c — u\\, which 
implies that p € dB c . Here again, by choosing 5 sufficiently small, we can make sure that the radius 
of curvature of the bump is at least ^, which means that the reach of the deformed hypersurface 
is still y = u- We can also make sure that the bump of p is disjoint from the bump of c since 
|| c — p\\ > 4=, and that the points of Lq \ {po} remain^ on M. It follows that B c is empty of points 
of L, where L is defined by L = Lq U {p} \ {po}- Since dB c contains u, v, w,p, B c is a Delaunay ball 
circumscribing [u,v,w,p]. Equivalently, c belongs to the Voronoi edge e dual to [u,v,w,p]. Note 
also that L is an (e — 5)-sparse (2e + <5)-sample of M. 

Since [u, v, w, p] is included in hyperplane z = 0, its dual Voronoi edge e is aligned with (0, 0, 1, 0), 
as illustrated in Figure [2] (right). This edge is incident to four Voronoi 2-faces, which are dual to the 
four facets of [u, v,w,p\. These 2-faces can be seen as extrusions, into the z-dimension (0,0, 1,0), 
of the edges of the Voronoi diagram of {u,v,w,p} inside hyperplane z = 0. Among these Voronoi 
edges, two lie above the plane t = A + f , and two lie below. As a result, in M 4 , two Voronoi 
2-faces incident to e lie above hyperplane t = A + |. These two Voronoi 2-faces do not intersect 
M, except at c and possibly at the bump of p. Now, the circumradii of the facets of [u,v,w,p] are 
at most ||c — u|| = ^]y^ 2 < rch(M), thus, inside hyperplane z = 0, Amenta and Bern's normal 
lemma [U Lemma 7] states that the edges of the Voronoi diagram of {u, v,w,p} make angles of at 
most arcsin < ^ with vector (0,0,0, 1). As a consequence, any Voronoi 2-face / incident to e 
in ]R 4 makes an angle of at most ^ with the plane passing through c, of directions (0, 0, 1, 0) and 
(0, 0, 0, 1) (note that the affine hull aff (/) intersects this plane along the line aff (e)). Since p lies 

away from this plane and only | above c, for sufficiently small 6 the Voronoi 2-faces incident to e 
lying above hyperplane t = A + 1 do not intersect the bump of p. As a consequence, they intersect 
M only at c, and therefore their dual Delaunay triangles are incident to exactly one tetrahedron 
of T> M (L), namely [u,v,w,p\. Hence, T> M (L) is not a closed hypersurface, and for this reason it 
cannot be homeomorphic to M. □ 

Observe that the example given in the proof corresponds to a degenerate case, since the Voronoi 
edge e dual to tetrahedron [u,v,w,p] intersects M tangentially at c. This degeneracy can be 
removed by inflating the bump of c infinitesimally, so that it intersects e twice and transversally, 
but still no other Voronoi edge. 

Note also that tetrahedron [u,v,w,p] is a sliver, since vertex p lies close to the affine hull of 
[u, v,w]. The original counter-example of [6] was designed to highlight the fact that the normals 
of slivers in the restricted Delaunay triangulation may differ significantly from the normals of the 
underlying manifold. This is not true for non-sliver simplices, as shown in Lemma 15 of [6]. Thus, 
the fact that [u,v,w,p] is a sliver in our counter-example is crucial. 

Theorem 2.2 For any positive constant fi < |, there exist a compact closed hypersurface M in M 4 
and an Q(e)-sparse 0{e)-sample L of M, with e = /x rch(M), such that D M (L) is not homotopy 
equivalent to M . The constants hidden in the f2 and O notations are absolute and do not depend 
on [i. 

Proof. Let A = — , and let S > be an arbitrarily small parameter. We begin our analysis 
with the example built in the proof of Theorem 12.11 We will modify M and L in such a way that 

1 They lie at least e away from po, and hence at least e — 8 away from (0, 0, 0, A). 
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Figure 3: Top: pentahedron \p,u,v,w, q] and the duals of \p,u,v,w], \p,u,v,q], and \p,v,w,q]. 
Bottom: the bump of c puvq and c pvwq (the bump of c has been slightly deflated). 

tetrahedron [p, u, v, w] will no longer belong to T> M (L) while its four facets will still. This will 
prevent T> M (L) from being homotopy equivalent to M. 

Consider point q = ( 1+ 2 V ^ , |,<5 2 , A + 5). The distance of q to hyperplane f = A is 5, which 
is arbitrarily small. Hence, as explained in the proof of Theorem 12.11 it is possible to deform M 
slightly by creating a small bump of radius of curvature at least A a t point ( l+ if^ , |, S 2 , A), so that 
M now passes through q while its reach remains y = jj. Moreover, since g lies farther than | from 
{p, it, to}, we can assume without loss of generality that its bump does not affect the positions of 
p, u, v, w. 

The circumcenter of pentahedron [p, u, v, w, q] is d = (i, i, 4-, A+ |), and its circumradius r' is 
less than e = 1 (for sufficiently small <5). It follows that every point of M lying in the ball B(c',r') 
is at distance 0(e) of {p, u, v, w, q}. In addition, q is farther than | from {p, u, v,w}. Therefore, if 
we modify L by inserting q and deleting all the points that lie in the interior of B(d , r'), L remains 
an 0(e)-sparse 0(e)-sample of M. Moreover, [p, u, v, w, q] is now a Delaunay pentahedron, whose 
dual Voronoi vertex is d . 

Note that line (c, d) is the affine hull of the Voronoi edge e dual to [p, u, v, w], and that d is an 
endpoint of e — see Figure [3] (top). Recall that, among the four 2-faces incident to e, two lie above 
hyperplane t = A+|. Let f puv and f pvw denote these two 2-faces. They are dual to triangles [p, u, v] 
and [p, v,w] respectively, since p lies above hyperplane t = A, which contains [n, v,w]. Moreover, 
f puv and f pvw are convex polygons, whose boundaries are two cycles of Voronoi edges that intersect 
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each other along e. In the cycle of df puv , one edge adjacent to e, noted e puvq , is dual to tetrahedron 
[p, u, v, q]. Similarly, in the cycle of df pvw , one edge adjacent to e, noted e pvwq , is dual to [p, v, u;, </]. 
Note that c' is an endpoint of both e puvq and e pvwq . Moreover, it can be easily checked that the 

line aft (e^q) passes also through point c puvq = (jj + - ; |j ^ + I + ^l+V^ )' wmle tne 
line aff(e pi , log ) passes through c pvwq = (5, 5 + <5 2 (<5 2 + 1), —5, A + f + 5(^ 2 + 1))- This implies that 
epuvq and ej, TO? make angles of 0(5) with hyperplane t = A + |. So, we are in a situation where 
tetrahedron [p, u, v, w] has a horizontal dual edge, while two of its adjacent tetrahedra, namely 
[p, u, v,q] and \p,v,w,q], have almost horizontal dual edges, as illustrated at the top of Figure El 

Since ||cp UV g p|| = ||cp Ui; g u\\ = \\c puvq — v\\ = \\c puvq q 1 1 < ||cp UU g ^11 j which is less than 
e = 1 for sufficiently small 5, we can modifjH L such that the ball B(c puvq , \\c puvq — q\\) contains 
no point of L in its interior, while L still remains an J7(e)-sparse 0(e)-sample of M. Similarly, we 
can assume without loss of generality that B(c pvwq , \\c pvwq — q\\) is a Delaunay ball. It follows that 
Cpuvq €E e P uvq and c pvwq £ Cp VW q- Since Cpuvq and c pvwq lie 0(5) away from each other, 0(6) above 
hyperplane t = A, and 0(e) away from L, we can deform M by creating a bump passing through 
Cpuvq and c pvwq , of height 0(<5) and radius of curvature at least while maintaining the points of 
L on M — see Figure [3] (bottom). Moreover, since c puV q and c pvW q also lie 0(e) away from c, we 
can assume without loss of generality that their bump does not touch the bump of c. It follows that 
tetrahedra [p, u, v,w], [p,u,v,q], and \p,v,w,q] belong to T> M (L), while M is still tangent at c to 
the Voronoi edge e dual to tetrahedron [p, u, v, w]. We call M + the current version of hypersurface 
M, and T> M+ (L) the Delaunay triangulation of L restricted to M + . 

Our last operation consists in deflating slightly the bump of c, such that M no longer intersects 
e, and thus [p, u, v,w] no longer belongs to V M (L), as illustrated at the bottom of Figure El Note 
however that \p,u,v,q] and [p, v,w,q] (and thus triangles \p,u,v] and [p, v,w]) are still in V M (L), 
since the bump of c pui)g and c pvwq is disjoint from the bump of c. Recall also that the Voronoi 
2-faces dual to [p, u, w] and [u,v,w] lie below hyperplane t = A + |, and that they make angles 
of at most ^ with vector (0,0,0, 1). Since the deflation of the bump of c is arbitrarily small, the 
Voronoi 2-faces dual to [p, u, w] and [u, v, w] still intersect M. As a consequence, the two triangles 
remain m V M (L). We call M~ the current version of hypersurface M , and T> M (L) the Delaunay 
triangulation of L restricted to M~. 

The result of these operations is that, although L is an 0(e)-sparse 0(e)-sample of both hyper- 
surfaces M + and M~ , whose homotopy types and reaches are the same as the ones of M, V M+ (L) 
and T> M (L) are different. Specifically, tetrahedron \p,u,v,w] is contained in T> M+ (L) but not in 
D (L), whereas its facets belong to both complexes. It follows that the Euler characteristics of 
T> RI+ (L) and D M (L) diffeJl, which implies that the complexes have different homotopy types. 
Therefore, at least one of them is not homotopy equivalent to the 3-sphere M. □ 

Witness complex. It is proved in [9] that the witness complex C w (L) may not contain V M (L) 
when the sets W, L are drawn from a smooth surface M such that W C M. However, we know from 
[3] that C W (L) is still included in V M (L) in this case. Below we prove that this latter statement 
no longer holds if M is a smooth manifold of dimension 3 or more: 

Theorem 2.3 For any positive constants /j,,u < |, there exist a compact closed hypersurface M 
in M 4 ; an Q.(e)-sparse 0(e)-sample L of M, and a 5-sample W of M, with e = \i rch(M) and 

2 For instance, we can simply delete the points of L that lie in the interior of the ball B(c pU vq, \\c puvq — q\\). 
Specifically, X (© M+ (£)) = x{V M ~ (L)) - 1. 
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5 = v rch(Af), such that C w (L) is not included in T> M (L). The constants hidden in the Q, and O 
notations are absolute and do not depend on \i nor v. Moreover, W can be made indifferently finite 
or infinite, and arbitrarily dense. 

Proof. Let M~, L, e, and c be defined as in the proof of Theorem 12.21 Recall that tetrahedron 
[p, u, v, w] does not belong to V M (L), whereas its facets do. We assume without loss of generality 
that c is not an endpoint of the Voronoi edge e, which means that the bounding sphere of the De- 
launay ball B(c, ||c — p\\) contains no point of L other than p, u, v, w. This condition can be ensured 
by an infinitesimal perturbation of the points of L \ {p, u, v, w}. Let d c = min p / gi \| p u „ w \ \\c — p'\\. 
This quantity is greater than ||c — p\\ since B(c, \\c — p\\) contains no point of L \ {p, u, v, w}. 

Consider any (finite or infinite) set of witnesses W C M~ such that, for each facet a of [p, u, v, w], 
W contains at least one point of M~ n V(u) (every such point witnesses a and its subsimplices). 
Assume further that W contains the top point of the bump of c (call this point c"). In the last 
stage of the perturbation of M described in the proof of Theorem 12.21 we slightly deflated the 
bump of c, such that c" lies strictly below c. Note that p is the vertex of \p, u, v, w] lying furthest 
away from c". Since the deflation is arbitrarily small, we can assume without loss of generality that 
||c-c"|| < \{d c - ||c-p||). This implies that the ball B(c", \\c" - p\\) C B(c, \\c - p\\ + 2||c - c"||) is 
included in the interior of B(c, d c ). As a result, B(c" , \\c" —p\\) contains no point of L \ {p, u, v, w}. 
Since p,u,v,w belong to B(c", \\c" — p\\), tetrahedron \p,u,v,w] is witnessed by c". And since the 
facets of \p, u, v, w] and their subsimplices are witnessed by points of W, \p, u, v, w] belongs to the 
witness complex C W (L). However, we saw in the proof of Theorem 12.21 that [p, u, v,w] does not 
belong to V M ~ (L). □ 

3 Conclusion 

We have proved that the structural properties of the restricted Delaunay triangualtion and witness 
complex on 1- and 2-manifolds do not hold on higher-dimensional manifolds. This implies in 
particular that the Delaunay-based approach to meshing and reconstruction is unlikely to work as 
is in higher dimensions. One possible way of getting rid of pathological cases is to use the sliver 
exudation technique of [5], which assigns weights to the vertices of the triangulation in order to 
remove slivers from the vicinity of the restrited Delaunay triangulation. This strategy has been 
successfully applied in [U [6] . 
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